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The influence of an external electric field on the current in charged systems is investi- 
gated. The results from the classical hierarchy of density matrices are compared with the 
results from the quantum kinetic theory. The kinetic theory yields a systematic treatment 
of the nonlinear current beyond linear response. To this end the dynamically screened 
and field-dependent Lenard-Balescu equation is integrated analytically and the nonlinear 
relaxation field is calculated. The classical linear response result known as Debye - On- 
sager relaxation effect is only obtained if asymmetric screening is assumed. Considering 
the kinetic equation of one specie the other species have to be screened dynamically while 
the screening with the same specie itself has to be performed statically. Different other 
approximations are discussed and compared. 



I. INTRODUCTION 

High field transport has become a topic of cur- 
rent interest in various fields of physics. In semi- 
conductors the nonlinear transport effects are ac- 
cessible due to femto - second laser pulses and 
shrink devices In plasma physics these field 
effects can be studied within such short pulse pe- 
riods One observable of interest is the cur- 
rent or the electrical conductivity which gives ac- 
cess to properties of dense nonideal plasmas || . In 
high energy physics the transport in strong electric 
fields is of interest due to pair creation . In order 
to describe these field effects one can start conve- 
niently from kinetic theory. Within this approach 
the crucial question is to derive appropriate kinetic 
equations which include field effects beyond linear 
response. 

At low strength of the external electric field one 
expects the linear response regime to be valid. 
Then the contribution of field effects to the con- 
ductivity can be condensed into the Debye- On- 
sager relaxation effect P-pd|] which was first de- 
rived within the theory of electrolytes ]T^lr| . De- 
bye has given a limiting law of electrical conduc- 
tivity |L2| which stated that the external electric 
field £ on a single charge Z — 1 is diminished in 
an electrolyte solution by the amount 



SE = E(1-—) 



(1) 



where e is the elementary charge, E the electric 
field strength, T is the temperature of the plasma 
and k is the inverse screening radius of the screen- 
ing cloud. This law is interpreted as a deceleration 
force which is caused by the deformed screening 
cloud surrounding the charge. Later it has been 



shown by Onsager 
corrected to 

SE = E (1 



1311 that this result has to be 
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if the dynamics of ions (Z — 1) is considered. 
While the linear response theory seems to repro- 
duce this Onsager result , the kinetic theory 
seems to support more the Debye result JTo| , p^7|JTl| ] . 
The correct treatment is a matter of ongoing de- 
bate. In this paper we will give the result beyond 
linear response for the statically and dynamically 
screened approximation. Here different approxi- 
mations of kinetic theory will be discussed and the 
one which leads to the closest form to the hydro- 
dynamical approximation (Onsager result) is pre- 
sented. 

The kinetic approach describes the time evolu- 
tion of the one particle distribution function within 
an external field E as 



d d 
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where the field dependent collision integral /[/, E] 
has to be provided by different approximations. In- 
tegrating this kinetic equation over the momentum 
k one obtains the balance of the current. For sim- 
plicity we assume that the distribution function 
can be parameterized by a displaced local equilib- 
rium one with a field and time dependent momen- 
tum /(k, t) = fo(k — p(E,<)) which is related to 
the current J as 
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if the charge is Ze, the density n and the mass 
TO. The balance equation for the field and time 



1 



dependent local momentum p(E, t) follows from 



dt P 
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)E = R(E)eZn3 (5) 



where the relaxation field SE(E) as well as the free 
conductivity R(E) follows from the field dependent 
collision integral. The total conductivity E = <rJ 
is then given by 



a(E) 
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The free conductivity R is the subject of intense 
investigations in the literature p). It is known that 
the Coulomb divergence for small wave vectors is 
omitted if screening is included and the divergence 
at large wave vectors is omitted by the De Broglie 
wavelength i.e. by the quantum effects. We will 
not consider the discussion of the free conductivity 
R here but concentrate on the relaxation field SE. 
The free conductivity can be obtained by the same 
considerations as will be outlined here. We want 
to point out that the relaxation field will turn out 
to be free of long wave divergences in the classical 
limit in contrast to the free conductivity R. 

First we recall the hydrodynamical approach 
starting from the classical Bogoliubov-Born- 
Green-Kirkwood-Yvon (BBGKY) hierarchy which 
results into an analytical formula for the classi- 
cal relaxation effect already reported [Is] . This 
result is then compared with the quantum ki- 
netic approach. We give a short rederivation of 
the field dependent kinetic equations in dynami- 
cal screened approximation from the Green's func- 
tion technique in Sec. III. Two approximations, 
the static screening as well as dynamical screening 
are presented. In the fourth section we will de- 
rive the field dependent current analytically. We 
present both the statically as well as dynamically 
screened treatment as analytical results. The clas- 
sical expressions for the statically screened result 
Jl7[ is compared with the classical result from hy- 
drodynamical approximation. The dynamical re- 
sult is then derived analytically too and it will be 
shown that only for asymmetric screening the hy- 
drodynamical result can be approached. In Sec. V 
we shortly discuss the physical limitation of field 
strengths for the local equilibrium assumption and 
the gradient approximation. Sec. VI summarizes 
and the appendix gives the calculation of some in- 
volved integrals appearing during the integration 
of the Lenard-Balescu equation. 



II. APPROACH BY CLASSICAL 
BBGKY-HIERARCHY 

The starting point for the classical considera- 
tions is the BBGKY hierarchy jn||o| which reads 
for the one - particle distribution function F a 

dF a dF e a -dF a 
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and the two - particle distribution function F a b 
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dt 
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with the external field E. S a describes a collision 
integral with some background which we will spec- 
ify later. This hierarchy is truncated approximat- 
ing that {H 

Fab = F a F b + g a b 

E a bc = F a F b F c + F a gbc + F b g ac + F c g a b (9) 

where <? a b(r a , i"b, v a , Vb) is the two-particle corre- 
lation function. 

Within the local equilibrium approximation 
we suppose a stationary (for example a local 
Maxwellian) distribution for the velocities in the 
one and two-particle distribution functions 



3/2 



fa(r,v,t) = n a (r,t) (^) 
g a b(r, r', v, v', t) = F ab - F a F b 
= h ab {r,r ,t) ^ \ 
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Here we have introduced the local one-particle den- 
sity and the local average velocity 



n a {r,t) = J dvF a (r,v,t) 

u a = — / dvv F a (r,v,t) 
n a J 



(11) 
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as well as the pair correlation function and the av- 
erage pair velocity 

h a b(r,r',t) = J dvdv'g ab (r,r',v,v',t) , 

w ab (r, r', t) = — [ dv dv'vg ab (r, r', v, v', t). 

n-ab J 

(12) 

Further on, we suppose that the particles interact 
with some background (e.g. neutrals or electrolyte 
solvent) by the collision integrals S a with the fol- 
lowing properties 



dvS a f a = 



dwS a f a 



bam a 



/ dvvS a g ab (r, r', v, v', t) = 1 fo afc w ab (13) 
J o a m a 

where b a is the mobility of particle of type a. This 
friction with a background serves here to couple 
the two - particle equations and will be considered 
infinitesimal small in the end. However, as we will 
demonstrate this yields to a symmetry breaking in 
the system which leads basically to different results 
than neglecting this friction. 

Fourier transform of the resulting two equations 
@ into momentum space and assuming a homo- 
geneous density n(r) = n we arrive at the coupled 
equation system 

ea-s / u u a v-^ 47r n b e a e b 

-_E.(v B -u.)/ = — + E^— 

b 

dk ik.(v a - w ab ) - . , > 

X / t^^t — / a (v a - W ab + U a )/l ab (k) 
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-e b k.(v b - u b ) / b ff oc (k)] + 5 a .g afc + S b g ab (15) 
with 

a i 



E = E - ^2 n b e b / dr b dv b 
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By multiplying the above equation system by 
1, v a , v b and integrating over the velocities we ob- 
tain the Onsager equation [fill 



Th ab (k)(l+i— a) + e Q $ b (-k) 



= -h 

with 



Th ab (k)(l - i^-a) + e 6 $ a (k) 



(17) 



fc 2 $ a (k) = 4vre a + ^ n c e c h ac (k) 

C 

fc 2 $ a (-k) = 47re a + ^ n c e c /i ca (k) (18) 

c 

for the two -particle correlation function h ab . Here 
we use 
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Let us already remark here that the friction with 
a background described by the mobilities b couple 
the two sides of the equation (0). If we had not 
considered this friction, Si = 0, we would have 
obtained that the left and the right hand side of 
( |l7| ) vanish separately. This will lead essentially 
to a different result even for infinite small friction. 
There is no continuous transition between these 
two extreme cases pointing to a symmetry breaking 
in the two treatments. Let us first discuss the case 
with background friction. 



A. With background friction 

The system ([l7]) for electrons, e e = e, and ions, 
e,i = —Ze, with charge Z reads expanded 

$ e (-k)+$ e (k) 



Th P 



Th e 



Th ie . = 



Thu = Ze 



^(-k)-Z|i$ e (k) 
l l + h- + l a{l+ b t Z) 
$i(k) - Z h t $ e {-k) 

$,(-k) + ^(k) 



(20) 



This we can solve together with (jig). First we 
calculate the effective field strength at the position 
of the electron in linear response the Onsager result 

Pi 



6E E 1 

1T E = -' 
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with k 2 =K 2 e {l + Z)= n ^ e n *> and 



_ 4Tv(e 2 n c + Z 2 e 2 ni) 
b e + Zbi 



(l + Z)(b e + b % )- 



(22) 



For single charged ions Z = 1 the influence of the 
mobilities drop out and we recover the result @). 

Since this result is independent of the mobilities 
one could conclude that this is an universal limit- 
ing law. However we will express two doubts here. 
As one sees for charges Z > 1 the result (|J) is 
approached only in the limit where the ion mobil- 
ities are much smaller than the electron mobilities 
bi/b e — > 0. This means of course that the electrons 
have different friction with a thought background 
than the ions. In other words there is an explicit 
symmetry breaking mechanism included by assum- 
ing such collision integrals with the background. 
Therefore we will obtain another solution if we con- 
sider no friction. 

The second remark concerns the limit of one- 
component system which one can obtain by setting 
Z = — 1. The Onsager result or hydrodynamical 
result with friction (|2l| ) leads to twice the Debye 
result (|l|) in this case but with opposite sign. Op- 
positely we will see in the following that the per- 
fectly symmetric treatment of the species without 
friction with a background will lead to a vanishing 
one component limit as it should. This again un- 
derlies the symmetry breaking if one assumes an 
infinitesimal small friction with a background. 

For completeness, we want to recall the expres- 
sion of the nonlinear Onsager result which 
is obtained from the limit bi/b e — > of the system 

r/tee + e ^(- k ) + ^ k ) =0 

h ei (T + ie^) + eifi(-k) = o(-j^) = 
k b e 

-h ie {T - i e^) - e<fi(k) = o(y ) = 



One obtains lP.4D.8l the result for Z = 1 
5E = 



(23) 
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with 
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The numerical values of this result will be dis- 
cussed in chapter ( IV C] ). 



B. Without background 

Now we reconsider the steps from ( |l5| ) to ( fl7| ) 
without friction with the background. We obtain 
that both sides of (O) vanish separately 



Th ab (k)(l + i^a) 

e 



Th ab (k)(l 



.e-b . 
i — a) 
e 



■e * t (-k) = 
e b $ a (k) =0. 



(26) 



Both equations have identical solutions h ab which 
can be easily verified using the symmetry h a b(k) = 
hba(—k). Together with ( p^| ) we can solve for $ e 
and the relaxation field is obtained instead of (p3) 



SE = -^±l(Z + l)BF N A (27) 

which takes for Z = 1 
SE e 2 K, 




(28) 



with 
F N (a) 
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^JA+{1 + Z)a 2 



■log- 
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(29) 



We see that the linear response result for Z = 1 
is twice the Debye result (0). For equal charged 
system Z = — 1 which would coincide with a one 
component plasma no relaxation effect appears as 
one would expect. In other words in a perfectly 
symmetric mathematical two - component plasma 
there is another relaxation effect than in a system 
which distinguishes the components by a different 
treatment of friction. The Onsager result (|24|) does 
not vanish for the limit of one-component plasma 
Z = — 1. This is due to the different treatment of 
ions and electrons there which assumes explicitly 
a two component plasma. Therefore the limit Z = 
— 1 does not work there. 

This result is quite astonishing. One would ex- 
pect that the limiting procedure which transforms 
the system (|l8|) into (|2^) would also lead to a 
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smooth transitions of the end results. However 
this is not the case. While the separate limit of 
b e ,i — > oo of (18) leads to ( |26| ) there is no possi- 
bility to transform the result ( pl| ) into the linear 
response result of (p8|). This underlines that due to 
even infinitesimal small friction assumed in obtain- 
ing ( ^l|) there occurs a symmetry breaking in the 
sense that the electrons and ions are not anymore 
symmetrically treated. 

This lesson we have to keep in mind when 
we now advance and investigate the systematic 
treatment by quantum kinetic theory. There we 
will find also complete different results when we 
use asymmetric screening compared to symmetric 
screening. Of course, we will not assume any phe- 
nomenological friction since the kinetic theory pro- 
vides for a systematic description of all occurring 
processes. Here we want only to point out that 
the above symmetry breaking is the main reason 
for the confusion in literature. Following the lin- 
ear response formalism an asymmetric treatment 
of two - particle correlation functions is used in 
that the electrons are statically screened . This 
seemingly innocent usage leads there to an occa- 
sional agreement for Z — 1 with the Onsager result 

Another advantage of the kinetic theory we want 
to point out here. The classical local equilibrium 
or hydrodynamical approximation does not lead to 
a mass dependence of the relaxation effect. This 
will be provided by the kinetic theory. 



III. QUANTUM KINETIC THEORY 

We will formulate the kinetic theory within 
gauge invariant functions not missing field effects. 
The most promising theoretical tool is the Green 
function technique Jl|,pl|-[23| . The resulting equa- 
tions show some typical deviations from the ordi- 
nary Boltzmann equation: (i) A collision broad- 
ening which consists in a smearing out of the ele- 
mentary energy conservation of scattering. This is 
necessary to ensure global energy conservation 24 1. 
(ii) The intra-collisional field effect, which gives ad- 
ditional retardation effects in the momentum of the 
distribution functions. This comes mainly from the 
gauge invariance. 

One of the most important questions is the range 
of applicability of these kinetic equations. Up to 
which field strengths are such modifications impor- 
tant and appropriate described within one-particle 
equations? In [g5| this question has been investi- 
gated for semiconductor transport. It was found 
that for high external fields the intra-collisional 
field effect becomes negligible. This range is given 



by a characteristic time scale of field effects rjj. = 
m?i/(eE-q) which has to be compared with the 
inverse collision frequency. This criterion is a pure 
quantum one. It remains the question whether 
there are also criteria in the classical limit. For a 
plasma system we will discuss in Sec. V that there 
is indeed a critical value of the field strength which 
can be given by classical considerations. 



A. Definitions 

In order to describe correlations in highly 
nonequilibrium situations, we define various cor- 
relation functions by different products of creation 
and annihilation operators 



G>(1,2) = < *(l)*+(2) > 
G<(1,2) = < #+(2)*(l) > 



(30) 



Here <> is the average value with the unknown 
statistical nonequilibrium operator p and 1 denotes 
the cumulative variables (ri, Si, t\...) of space, 
spin, time etc. The equation of motion for the 
correlation functions are give n in the form of the 
Kadanoff-Baym equation P^ , f27j . p3[ 



i(G ( 7 1 G < - G < G(7 1 ) = i (G^ST 
-i (Y, R G < - G<^ A ) 



<nA\ 



E < G 



(31) 



where the retarded and advanced functions are in- 
troduced as A R (1, 2) = -iQ(t 1 -t 2 )[A > ±A<] and 
A A (1,2) = iQ(t 2 -t 1 )[A > ±A < }. Here operator no- 
tation is employed where products are understood 
as integrations over intermediate variables (time 
and space) and the upper/lower sign stands for 
Fermions/Bosons respectively. The Hartree- Fock 
drift term reads 



d 



2 

xi 



S ffF (ll')J a(i-i') 
(32) 

with the Hartree Fock self energy 

= ^(ri-ri) j rfr 2 U(r 1 -r 2 )G < (r 2 t' 1 r 2 < 1 ) 
+V(ri - r' 1 )G < (r 1 i 1 r' 1 <' 1 )) 6(fi - t' x ) (33) 
where G(r2,ti,r 2 , t\) = n(r 2 ,ii) is the density. 

B. Gauge invariance 

In order to get an unambiguous way of con- 
structing approximations we have to formulate our 
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theory in gauge invariant way. This can be done 
following a procedure known from field theory (28) . 
This method has been applied to high field prob- 
lems in . With the help of the Fourier transform 
of an arbitrary function G(x,X) over the relative 
coordinates x = {r% — ri, t% — 1\) = (r, r) with the 
center of mass coordinates X = ((r2 + ri)/2, (<2 + 
t\)/2) = (R, t) one can introduce a gauge-invariant 
Fourier-transform of the difference coordinates x 

G{k,X) = { dxG(xX) 



xexp^ / d\x )i [k tI + -A fl (X + Xx)}} . (34) 



For constant electric fields, which will be of interest 
in the following, one obtains a generalized Fourier- 
transform 



G(k,X) = / dxeK [x " 



G(x,X), 



where the \ function was chosen in such a way that 
the scalar potential is zero = (0, — cEt). There- 
fore, we have the following rule in formulating the 
kinetic theory gauge-invariantly 

1. Fourier transformation of the 4-dimensional 
difference- variable x to canonical momentum 
P- 

2. Shifting the momentum to kinematic mo- 
mentum according to p = k — eEt. 

3. The gauge invariant functions G are given by 

G(p, t) = G(k - eEt, t) = G(k, t) 

= G(p + eEt,t). (35) 



Here /(p,t) = G < (p,R, t, r = 0) denotes the 
Wigner distribution function and we suppress the 
center of mass coordinates. {,} is the anti- 
commutator understood that the r variable at 
the first place comes with a minus sign respec- 
tively. This equation is exact in time, but accord- 
ing to the assumed slowly varying space depen- 
dence we have used gradient expansion for space 
variables and dropped all R-dependence for sim- 
plicity. This criterion is discussed in the last sec- 
tion (96). With the help of the gauge invariant for- 



mulation of Green's function ( pl^ ), we can write 
the kinetic equation ( ^6| ) finally in the following 
gauge-invariant form 

t-t 

^/(k,t)+eEV k /(k,t)= J dr 



G > (k- ^r, r, t~), £< (k- ^r, -r, t- J 
v 2 ' ' 2' v 2 2 



-^<(k-^r,r,t-I),E>(k-^r,- 



2' 



(37) 



This kinetic equation is exact in time convolutions. 
This is necessary because gradient expansions in 
time are connected with linearization in electric 
fields and consequently fail f32fl . The gradient ap- 
proximation in space has been applied assuming 
slow varying processes in space. This corresponds 
to the limit of a weakly coupled plasma, which we 
employed already in Section [ll]. Please remind that 
due to Coulomb gauge we do not have space inho- 
mogeneity by the electric field. 



We shall make use of these rules in the following 
sections. In Q this procedure has been general- 
ized for two - particle Greens functions and leads to 
the field - dependent Bethe - Salpeter - equation. 



C. Equation for Wigner distribution 

In the relative and center of mass coordinates 
the time diagonal part of ([|l|) reads jn| 



d_ 

dt 



f(p, t) 



L — LQ 

J dT[{G > (p,<-^,-r),E<(p,<-^,r)} 
o 

- {G<(p,i-1 - r ),£>(p,i-l r)}' 



(36) 



D. Spectral function 

The spectral properties of the system are de- 
scribed by the Dyson equation for the retarded 
Green function. A free particle in a uniform elec- 
tric field, where the field is represented by a vector 
potential E(i) = — i A (t) leads to the following 
equation 



* 9 f 
lh dt - £(P 



■-A(t)) 



G*(p,tt , ) = 6(t-t'). 



(38) 



This equation is easily integrated p3lt 



G£(p,tt') = -iQ(t-if)exp 



Jdue(p — -A(u)) 

t 

(39) 
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For free particles and parabolic dispersions, the 
gauge invariant spectral function |33|]34[| follows 



G<(p,t,T) = f(p,t 



-)A(p,r,i). 



(42) 



Ao(k,w) 



n f , ( k 

2 / arcos lot 
o 

2tt , ( k 2 /2m - Tilo 
— Ai — - 



e 2 E 2 



2mh 24mh 



(40) 



where Ai(x) is the Airy function ]|5j and = 
(h 2 e 2 E 2 /8m) 1 / 3 . It is instructive to verify that 
( f40| ) satisfies the frequency sum rule J <LjA (uj) = 
2ir. The interaction-free but field-dependent re- 
tarded Green's function G^ can be obtained from 
the interaction-free and field-free Green's func- 
tion by a simple Airy transformation |36|. This 
is an expression of the fact that the solutions of 
the Schrodinger equation with constant electric 
field are Airy- functions. The retarded functions 
can therefore be diagonalized within those eigen- 
solutions [B7LE9[. It can be shown that (Eo|) re- 



mains valid even within a quasiparticle picture 36 1 , 
where we have to replace simply the free dispersion 
k 2 1 2m by the quasiparticle energy ej,. 



E. The Problem of the ansatz 

In order to close the kinetic equation ([36]), it 
is necessary to know the relation between G > and 
G < . This problem is known as an ansatz and must 
be constructed consistently with the required ap- 
proximation of self-energy. Assuming the conven- 
tional KB ansatz [^6| we have a relation between 
the two time Green functions and the distribution 
function 

G<(k,w,r,t) = yl(k,w,r,t) /(k,r,t) 
G>(k,w,r,t) = yl(k,w,r,t) (1 T /(k,r,t)). (41) 

This is quite good as long as the quasi-particle pic- 
ture holds and no memory effects play any role. As 
we shall see, the formulation of kinetic equations 
with high fields is basically connected with a care- 
ful formulation of retardation times. Therefore, 
the simple ansatz, called KB ansatz fails. 

Another obscure discrepancy is the fact that 
with the old ansatz, one has some minor differences 
in the resulting collision integrals compared with 
the results from the density operator technique. 
With the old ansatz, one gets just one half of all 
retardation times in the various time arguments 
pl|p3[ . This annoying discrepancy remained ob- 
scure until the work of Lipavsky, et al. |3j| where 
an expression is given for the G < function in terms 
of expansion after various times. We can write in 
Wigner coordinates 



This generalized- Kadanoff- Baym (GKB) - ansatz 
is an exact relation as long as the selfenergy 
is taken in Hartree- Fock approximation. To- 
gether with the requirement of gauge invariance of 
Sec. |1II B| and using the quasiparticle spectral func- 
tion Q40D with quasiparticle energies instead of 
k 2 /2m, the GKB ansatz finally reads 

e 2 E 2 



G < (k,r,r,t) = exp < -- e k r 



2 2 



24m 



(43) 



In order to get more physical insight into this 
ansatz one transforms into the frequency represen- 
tation 



pEt t 
G<(k,w,r,t) = 2 / dr/(k -—,«--) 



xcos lot — e(k, r, £) — 



t e 2 E 2 



% 2AmU 



(44) 



Neglecting the retardation in / one recovers the or- 
dinary ansatz ( 41 ) with the spectral function (|io| ) . 
The generalized ansatz takes into account history 
by an additional memory. This ansatz is superior 
to the Kadanoff-Baym ansatz in the case of high 
external fields in several respects p9| : (i) it has 
the correct spectral properties, (ii) it is gauge in- 
variant, (iii) it preserves causality, (iv) the quan- 
tum kinetic equations derived with Eq.(^2|) co- 
incide with those obtained with the density ma- 
trix technique [[iol , [i"l|]3~i|l , and (v) it reproduces the 
Debye-Onsager relaxation effect . 

Other choices of ansatz can be appropriate for 
other physical situations. For a more detailed dis- 
cussion see (42]. 



F. Kinetic equation in dynamically screened 
approximat ion 

For Coulomb interaction it is unavoidable to 
consider screening if one does not want to obtain 
long range or short wave vector divergences. To 
obtain an explicit form for the kinetic equation we 
have to determine the selfenergy 'S >,< . The dy- 
namically screened approximation is given by ex- 
pressing the self energy by a sum of all ring dia- 
grams. The resulting kinetic equation is the quan- 
tum Lenard - Balescu equation, which has been 
derived for high fields in ]3(| . We give this approx- 
imation in exact time convolutions. The selfenergy 
is given in terms of the dynamical potential V 
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E<(M,f) = / ^3V<(q,*,f)G<(k-q,t,f) 

(45) 

where the dynamical potential is expressed within 
Coulomb potentials V a b(q) 



V< (q, t, = £ ^(q)r< (q, t, i')Ka(q) (46) 

via the density-density fluctuation 
£< b (q,t,t') = <5 o6 J didt 

x (f)- 1 (q,*,i)L< (q,t,t) (^r 1 (q,?,f'). (47) 
Here L is the free density fluctuation 

J&(q,*,f) = 1 ^^G<(p,M')G>(p-q,t',i) 

(48) 

and £ r/,<1 the retarded/ advanced dielectric func- 
tion 

r/«(q, i, i') = *(t - *') ± i0[±(t - 0] ^ V bb (q) 

b 

x(i>(q,i,t')-i < (q,i,t'))- (49) 

One easily convince oneself that this set of equa- 
tions is gauge invariant. 

We can directly introduce this set of equations 
into the equation for the Wigner function ( |37| ) and 
obtain after some algebra for the in-scattering part 
of the collision integral 

oo 

b 

e a Eqr 2 ~ 



X cos 



(4- 



•i " ~" 2m, 
x/ (k-q-e Er,t-T)(l-/ (k-e Er,t-r)) 
L^(q,w,t- |r) 



|£(q,w,t- §t)| 
with the free density fluctuation ( fls] 

oo 

dp 



(50) 



X cos 



(w - 4 + 4+ g ) T + 



fcEqr 2 

2m b 



x/ b (p + q, t - lr)(l - / b (p, i - ir)). (51) 

The out-scattering term 7 out is given by / «-» 1 — /. 
Here we used the ansatz (O) and have employed 



the approximation t ± It sa i in the density fluc- 
tuation (47) which corresponds to a gradient ap- 
proximation in times for the density fluctuations. 
Since the center of mass time dependence is car- 
ried only by the distribution functions in (|47| ) , this 
approximation is exact in the quasistationary case 
which we investigate in the next section. All inter- 
nal time integrations remain exact. Of course, for 
time dependent phenomenae we have to question 
this ap pro ximation. 

Eq. ( |50| ) represents the field dependent Lenard- 
Balescu kinetic equation [ 36| which was here 
slightly rewritten and which form will turn out 
to be very convenient for the later analytical in- 
tegration. Other standard approximations like the 
T-matrix J3(| approximation resulting into a field 
dependent Bethe-Salpeter equation can be given. 



1. Kinetic equation in statically screened 
approximation 

Using the static approximation for the dielec- 
tric function £ (q, 0, t) in (|50|) ; the kinetic equation 
for statically screened Coulomb potentials in high 
electric fields appears 1 10 3l|j36] ] 



d d 

df fa + eE dK fa = ^ Iab 

b 

2(2s 6 + l) f dKdk h dk' h 

lab = J (2-KTlf ( k a + K b -k a -k b ) 

X {fa' Ml ~ fa)(l ~ fb) - fafb(l ~ /„')(! " fb>)} 

OO 

xV s 2 (k a -k'J J drcos{(e Q + e 6 -e' a -e b )J 



Er 2 



2h V m a 



e b k b 

m b 



e bK 



(52) 



with fb = fb(kb — ebEr,T — r). The potential is 
the static Debye one 



V s (p) = 



4:Tre a e b h 



p 2 + h 2 n 2 

and the static screening length k is given by 

2 _ 4 7 re 2 .n e 
hi — y 



(53) 



(54) 



in the equilibrium and nondegenerated limit. Here 
T c is the temperature of specie c, charge e c , spin 
s c and mass m c respectively. 

If we had used the conventional Kadanoff and 



Baym ansatz (41 ) we would have obtained a factor 
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1/2 in different retardations Jul] . This would lead 
to no relaxation effect at all |10|| . Furthermore it is 
assumed, that no charge or mass transfer will occur 
during the collision. Otherwise one would obtain 
an additional term in the cos - function propor- 
tional to r 3 . 

Two modifications of the usual Boltzmann colli- 



sion integral can be deduced from ( 52 ) : (i) A broad- 
ening of the ^-distribution function of the energy 
conservation and an additional retardation in the 
center-of-mass times of the distribution functions. 
This is known as collisional broadening and is a 
result of the finite collision duration ■ This ef- 
fect can be observed even if no external field is 
applied. It is interesting to remark that this colli- 
sional broadening ensures the conservation of the 
total energy [p3. If this effect is neglected one 
obtains the Boltzmann equation for the field free 
case, (ii) The electric field modifies the broadened 
6- distribution function considerably by a term 
proportional to t 2 . This broadening vanishes for 
identical charge to mass ratios of colliding parti- 
cles. At the same time the momentum of the dis- 
tribution function becomes retarded by the elec- 
tric field. This effect is sometimes called intra- 
collisional-ficld effect. 



IV. FIELD EFFECTS ON CURRENT 

We are now interested in corrections to the parti- 
cle flux, and therefore obtain from (^2|) the balance 
equation for the momentum 

d 



dt 



< k a > 



-n a e, 



,E 



E < k » ia B > 



(55) 



Here we search for the relaxation field (|5|) 
which will be represented as renormalization of 
the external field E similar to the Debye-Onsager- 
Relaxation field in the theory of electrolyte trans- 
port JI^-[l^] . This effect is a result of the deforma- 
tion of the two-particle correlation function by an 
applied electric field. 

To proceed we assume some important restric- 
tions on the distribution functions. First, we as- 
sume a nondegenerate situation, such that the 
Pauli blocking effects can be neglected. Second, to 
calculate the current for a quasistationary plasma 
we choose Maxwellian distributions analog to din) 



fi(p) = 



2s, 



1 



A- exp 



P 



2m, ; Ti 



(56) 



with the thermal wave length A 2 = 2Trh 2 /(m l T l ), 
the spin Si and the partial temperature Ti for 
species i which can be quite different e.g. in a 
two - component system. 



A. Statically screened result 

Before we present the result for the dynamically 
screened approximation we want to give the static 
result. The momentum conservation in ( |52"| ) can 
be carried out and we get for the relaxation field 



■£(2s + l)(2s 6 +l) 



oo 



x cos 



kq qQ 



Eq 




2h 


m b m a 



2m a h m a h nibh 
x {/a(k)/ b (Q + | ) - / a (k + q)/ b (Q - f )} (57) 

where we have shifted the retardation into the dis- 
tribution functions. The second part of the distri- 
bution functions can be transformed into the first 
one by putting k + q — > k and q — > q with the 
result 

n a e a ^-'E = ^2 ^^Lo [dkdqdQ f b (Q) /„(k) 



, h 2 (2Trhf 



x V 2 (q) q J dr cos 



Eq 

~2h 



o 

mb m a 



q 2 _ _kq_ qQ \ 

2[ih m a Ti rribh J 



(58) 



with the reduced mass pT 1 = l/m a + I /mi,. The 
angular integrations can be carried out trivially 
and we get 

n a e a — E= -^/i 



h = 



e-b e a 
m b m a 



xsin(g)/ 2 H/ 2 [6] 



(59) 



withjs(x) = (x cos x — sin x)/x 2 . The two integrals 
over the distribution functions I2 can be done with 
the result 



h[a] 



fim a (2sg + 1) 
qr 



dkkf a (k) sin(-^-) 
m a h 



(60) 



and correspondingly I2 [b] . We now introduce the 
new variables 
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q = 2y^J fiT ab 
, 2T ab T 



1 



b = 7T 



-n 



AT 



3/2 



a b 



eb _ _ea_ 

m b m a 



(61) 



and obtain 



(62) 

Using the screened Debye potential ( ^3|) we finally 
obtain 

J l ~ 77; *3 



(z 2 + l) 2 





J7 • / ,2^ m ( Z ^C) -z 2 / 2 

dl]s{xzl ) — e 



(63) 



Therein we used y = z( and Z = i£ with the quan- 
tum parameter 



2 ^-2 



< 2 



AflT a b 

and the classical field parameter 

E ( m a m b 



x = — = 



(■b - 



( 2T a bK \m a + m b m a + mb 



(64) 



(65) 



With this form ( p3|) we have given an extremely 
useful representation because the field effects, con- 
tained in x, are separated from the quantum ef- 
fects, which are contained in £. The integral in 
( |63| ) can be performed analytically in the classical 
limit ( — > 0. For the more general quantum case 
with arbitrary £ the linear and cubic field effect 
can be given analytically and are discussed in 17 1. 
We will not discuss them here. 

Performing the classical limit £ — + one obtains 
from (pi that M 



he 

F(x) 



7VX 

24 
3 

~~ 2 



F(\x\) 
3 — x 



x 



1 + X X 



ln(l + x) 



(66) 



Introducing the classical result (|66|) into ( |63| ) we 
find from ( |59| ) and (55) the following relaxation 
field 

d ( SE \ 

— < k a > -n a e a E ( 1 + J = 72 a e a JR(E) 



with 

SE a 
E 



6k V 



4 2 ^k. _ a. 



iab / Tj, , T a 



(67) 



-F(\x\) (68) 



and x from (65). We see that for a plasma consist- 
ing of particles with equal charge to mass ratios, 
no relaxation field appears. The link to the known 
Debye- Onsager relaxation effect can be found if we 
assume that we have a plasma consisting of elec- 
trons (m e , e e = e) and ions with charge = eZ 
and temperatures T e = Ti = T. Then ( |68[ ) reduces 
to 

SE, = «e» Z(1 + ^Z) J eE Z (l + %Z) \ 
E 6T(1 + Z)(1 + ^) ^« l+sj ; 



6T 



forZ = 1. 



3kT 
2e_E 



(l/£) 2 



(69) 



This formula together with the general form ([38]) is 
the main result of this chapter. It gives the classi- 
cal relaxation effect for statically screened approx- 
imation up to any field strength and represents a 
result beyond linear response. We see that in the 
case of single charged heavy ions the Debye result 
(0) is underestimated by a factor of two. 



B. Dynamically screened result 

The calculation of the current with the collision 
integral for dynamically screened potentials ( [so] ) 
can be performed analytically as well. For the 
quasistationary condition we can calculate the fre- 
quency integral in ( |50| ) analytically using the iden- 
tity (44j for the classical limit o(h) 

fdu)H(uj) -I . , H(0) ( 



£(q,0) 

(70) 

where we set H{w) = w/Im£ and which relation is 
proven in appendix |X| We will employ only classi- 
cal screening. The quantum result for screening is 
more involved and not yet analytically integrable. 

Observing that for the dielectric function ( fi"9| ) 
together with (p6|) holds 
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lim — = tt 



'•'6 



(71) 



5E d * n 



ienK 



m a m b 
T a T b 



with the partial screening length = 4ire 2 n b /T b 
and the partial thermal velocity v 2 = 2T b /m b , we 
obtain for the current ( |59| ) after similar integra- 
tions as in chapter IV A instead of (pS3) 



3/2 , 



T, 



3/2 



+ o(E). (76) 



'mi 



7 dyn _ SK z eie(n a n b y/m a m b jdyn 



The difference to fl68| ) becomes more evident if we 
consider again only electrons and ions with equal 
temperature 



xdyn _ 



= / dz- 



1 o 



x - cos[M 6 C^ 2 + Bzl 2 x] cos[M a Chz 2 - Azl\x\. (72) 
Here we used the same dimensionless variables as 



8E^ n 
E 



.2 2Z(l 



6T 



(Z- 



-) 



o(E). (77) 



in chapter IVA and the quantum parameter (]6J). 



Further we abbreviated A = ^jr-, B 



2^T 



M b 



2fiT 



m a T a ' y m b T b ' 

We like to remark that we neglect any field de- 
pendence on the screening £ itself here. As pre- 
sented in ]39t a field dependent screening function 
can be derived. However, this field dependence 
gives rise to a field dependence starting quadrati- 
cally and will be not considered in this work. 

The classical limit of ( |7^ ) can be performed 
again by ( — ► 0. We obtain 



i: 



dyn 



1 



-AM a I[\A\,\B\] - (a~b) 



(73) 



with the remaining 3-dimcnsional integral 

oo 1 



I[A,B] 



dz- 



z 2 + l 



dx 



A 2 x 2 



die z 1 cos(Bzfx) 



(74) 



1. Linear response 

The linear response can be read off directly from 
7 J) and is given by I[0, 0] of (||). We obtain 



r dyn 



TT 3 / 2 

~Y2 



(M a A - M b B) 



(75) 



The differences to (|6^) are obvious in the different 
mass dependence. This result overestimates the 
Debye result by a factor of two. 



2. Complete classical result 

Now we are able to present a complete field de- 
pendence beyond linear response. The integral 
( [74] ) can be done analytically, which is sketched 
in appendix |b[ The result reads 



3 /2 

I[A,B} = ^—1[A,B] 



6 



4A(l -y/T+B A 2 +\og(l-A 2 ) 



B 



ArcTanh( 1 ) - ArcTanh(^±|) 



-ArcTanh( 



1 - 



) + ArcTanh( 



x/T+B • 



'A 



(78) 



We obtain for © 
5E^ n 4e a 7TK 



E 



^ k-2 rn a 



E 2 m a m b 
n b e b 



3/2 — 



T^3/2 



T a T b 
1[B, A] 



(79) 



Expanding ( JTq ) in powers of -E we recover (|76|). 
Once more we choose the case of electrons and ions 
with equal temperature and obtain 



and the linear relaxation field (Bq) takes the form 



5i? dyn _ Ke 2 2Z(1[A,B] + ^Z1[B,A}) 

E Qe T (z+,[™^) 



(80) 



11 



For single charged ions and big mass differences we 
can further simplify to 



<5£' dyn 



= -— T — 



9 

e n e 



T[x] 



6T kT 1 
2 + o(E) 

( 2 (-2a; + 3 (-1 + ^/T 



X 



V2 ( -ArcTanh(i) 
V v2 



ArcTanh( 



V2 



log(l - x 2 ) 
V2 



= 2 + o{E). 



(81) 



This result will be compare with the statically 
screened result (|69| ) and the hydrodynamical re- 
sult (pit) in section IV C . Here we remark already 



that the Debye result is twice overestimated here. 



Thermally averaged dynamically screened 
result 



We will now give an approximative treatment 
of the dynamical screening used in H. This ap- 
proximation consists into the replacement of the 
dynamical screening in the collision integral ( |50| ) 
which is £{oJ,q)~ 2 by (1 + K 2 V a Jq) / Att)' 1 . This 
represents a thermal averaging |44j of £~ 2 which 
can be proven easily with the help of a p pendi x W . 
We obtain the relaxation effect of ( |6q ) and ( |69f ) 
but with a different field function F 



F d ^(x) = -4 



2-x ln(l + x) 

x 

2 + o{x) 
l+o{l/x) 2 ' 



(82) 



Therefore the relaxation effect (|69|) in linear re- 
sponse for single charged ions takes the form of (]l|) 
and twice the static screened result ( |6^ ) and half 
of the dynamical screened result ( |8l| ) . 

As we see from figure [j] the different approxima- 
tions lead to very different results. The statically 
screened result (|6^ ) underestimates the Debye re- 
sult by factor 2 which is corrected by the ther- 
mal averaged treatment of the screening. If we 
calculate instead the complete dynamical screened 
result J77| ) or (|so| ) we obtain twice the Debye re- 
sult (Q) and the thermal averaged screened result. 
However there is a complete different charge de- 
pendence. We have to observe that the perfectly 



symmetric treatment of screening does not repro- 
duce the hydrodynamical result which is the On- 
sager result (0) for linear response. 



dyn. screened 

hydrodyn. with background 

statical screened 

therm, av. dyn. screened 

hydrodyn. without background 
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eE/KT 

FIG. 1. The nonlinear Debye-Onsager relaxation ef- 
fect vs. scaled electric field for an electron and single 
charged ion system. The hydrodynamical approxima- 
tion ( |24| ) leads to the Onsager result (Q) for small field 
strength 2 — \/2. The statically screened result ( p^ ) 
or ( |69[ ) leads to half the Debye result ([j]). The ther- 
mally averaged approximation of the dynamical screen- 
ing ( p2] ) leads to the Debye result while the full dy- 
namically screened approximation (^) leads to twice 
the Debye result. Also the hydrodynamical result ( pj| ) 
without background leads to twice the Debye result. 



D. Asymmetric dynamical screened result 

We want now to proceed and ask under which 
assumptions the Onsager result (||) might be re- 
produced. Following the results we saw from the 
hierarchy we have consequently to treat the elec- 
trons (specie a) and ions (all other species) asym- 
metrically. This we will perform in the same spirit 
as Onsager in that the ions have to be treated dy- 
namical (as before) but the electrons are screened 
statically. 

This means we consider as the potential not the 
bare Coulomb one but a statically screened Debye 
potential for specie a. The ions (all other species) 
will then form the dynamical screening. In com- 
parison with the chapter before we can perform 
all steps analogously except two modifications, eq. 
and 



we have now 



Urn 



First we observe that instead of ([71 



qv b (q 2 + h 2 K 2 e ) 



Im£ (q,ui) y/irh 3 Kl 
which leads to a replacement of the sum 



(83) 
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K 



(84) 



in the for-factor of ([72) and (|7^ ) . This leads in the 
limit of big mass differences to a for-factor in ( |77j ) 
and (80) respectively 



modification I 



Z 



(85) 



The second modification is that in (|t^) one has to 
replace 



1 + z 2 



1 + z 2 q + z 2 

~2 



with 



lg 



1 1 + z 2 



(86) 



(87) 



This shows that in the end results (|7^), Eq. (|8 
has to be changed 



modification II : 
1[A, B] ■ ^ 



-1[A,B}. 



A B 



Particularly we obtain for the linear response result 
( |77j ) where for electron-ion plasma q = 1/{Z + 1) 
and 



5E as y 6E & y n Zq 



E 



2 



HC 



Zq 



3T 



o(E) 



(89) 



which agrees with (|2l]) if we consider that the mo- 
bilities are very different bi/b e — > in (^2|). 

The same result we obtain from the thermally 
averaged result ( |66| ) since there appears no such 
function as ( |83] ) and therefore the modification I 
of (^5) does not apply but solely the modification 
II ofTfl). We therefore obtain @ but 



F a dyn (x) = 2F dyn (:r) - V2F dyn (V2x) 
2-V2 + o(x) 



3_ 

2.v 



o(l/ x) 



(90) 



with F dyn of (|8^) . The linear response leads then 
exactly to the same result as from the dynamical 
screening (|89"|), i.e. the Onsager result with the 
same charge dependence. 



The fact that we reproduce the classical Onsager 
result with the same charge dependence can be 
considered as very satisfactory. The more since 
we have seen how many different considerations 
are possible. Please note that the special case 
Z = 1 could lead occasionally to a seemingly agree- 
ment between different treatments. We think that 
the charge dependence incriminates different treat- 
ments. 



dyn. asymmetrical screened 

v. x hydrodyn. with background 

thermod. av. asymmetrical screened 




eE/KT 

FIG. 2. The different asymmetric screened approxi- 
mations for the relaxation effect versus field strength. 
The hydrodynamical result ( |24| ) is compared with the 
thermally averaged asymmetric screened result of ( po| ) 
as well as the asymmetric screened one of (fji|). 

In fi gur e ^ we see that the asymmetrical screened 
result ( p0| ) with ( |8S| ) approaches the hydrodynami- 
cal or Onsager result (^) rather well for small fields 
while it is too low at high fields. On the other hand 
the thermally averaged symmetrical screened re- 
sult (|9^) agrees with the hydrodynamical approx- 
imation (|2J) in the low and high field limit. Why 
the hydrodynamical result cannot be reproduced 
completely within the kinetic theory remains still 
a puzzle. Probably the remaining difference is due 
to the neglect of the field effect on the screening 
itself M. 



V. RANGE OF APPLICABILITY 

During the derivation of the quantum kinetic 
equations there has been assumed the gradient ap- 
proximations which restricted the spatial gradients 
of the system. Here we want to discuss up to which 
field strength this assumption is justified. 

The electric field is limited to values x « 1 for 
x from (^) . This can be deduced from the expres- 
sion for the dynamical screened result (|8l|). The 
expression has a remove-able singularity at x = 1. 
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Therefore we see a smooth curve. Nevertheless this 
is the field strength where something is happening. 
For equal masses and temperatures of plasma com- 
ponents this condition translates into 



E < 



kT 



(91) 



We interpret the occurrence of such singular point 
that no thermal distributions are pertained in the 
system. Then we have to take into account non- 
thermal field dependent distributions which have 
been employed to study nonlinear conductivity 



The condition (|9_y) allows for different physical 
interpretations. Within the picture of the screen- 
ing cloud we can rewrite ( |9l| ) into 



eE < m-^-. 



(92) 



This means that a particle moving on the radius of 
the screening cloud td = 1/ 'n with thermal velocity 
v% h = T/m should not be pulled away by the acting 
field force. We can discuss this limit also via the 
energy density which can be reached in a system 
by the applied field. We can reformulate once more 
the condition (Bfi) to find equivalently 



47T 



< nT. 



(93) 



This means that we have essentially non-thermal 
effects to be expected if the energy density of the 
field becomes comparable with the thermal energy 
density. 



The validity criterion (91) can now be used to 



check the weak space inhomogeneity which has 
been assumed during our calculation. Quasi- equi- 
librium in charged systems with external fields can 
only be assumed if the field current is accompanied 
by an equivalent diffusion current 



jfieid = e/j,En = - jdiff = eD 



dn 
dx 



(94) 



using the Einstein condition /i = eD /T one gets 



eE 



T- — 

n dx 



(95) 



Combining this elementary consideration with our 
condition ( |9l| ) we obtain a limitation for space gra- 
dients 



dn 



d (k x) 



< n 



(96) 



where our treatment of field effects and local equi- 
librium is applicable. 



VI. SUMMARY 

The nonlinear relaxation field of a charged sys- 
tem under the influence of high electric fields is in- 
vestigated. The local equilibrium or hydrodynam- 
ical approach starting from the classical BBGKY 
hierarchy is compared with the results from the 
quantum kinetic equations. 

We come to the same conclusion considering 
the hydrodynamical approximation or the kinetic 
theory that a perfectly symmetric two compo- 
nent plasma will lead to a different relaxation ef- 
fect than the case where we consider the mov- 
ing charge asymmetrically from the screening sur- 
rounding. In the hydrodynamic approach this has 
been achieved by friction with a background, in the 
kinetic approach we have realized it due to asym- 
metric screening. Within this asymmetric treat- 
ment the limit to a one component plasma which 
would be to set the ion charge to Z = — 1 leads to 
a non-vanishing finite quantity. Oppositely in the 
perfectly symmetrical treatment this limit vanishes 
in that the relaxation field vanishes as it should. 
The perfectly symmetrical treatment of species in 
the system leads to twice the Debye result different 
from the Onsager result in linear response. 

Different approximations of the kinetic ap- 
proaches are compared and discussed. We found 
that the symmetrical treatment of species as well 
as the asymmetrical treatment leads to the same 
corresponding results as the hydrodynamical ap- 
proach for linear response. But for higher field 
strengths there appear minor differences which are 
probably due to the neglect of the field dependent 
screening itself. The thermally averaged approx- 
imation of screening has the advantage to agree 
for low and high fields with the hydrodynamical of 
local equilibrium approach. 



ACKNOWLEDGMENTS 

The author wants to thank warmly Jens Ortner 
who has contributed the hydrodynamical approach 
with background. Also many comments regarding 
the manuscript are gratefully acknowledged. Last 
but not least Gerd Ropke is thanked for stimulat- 
ing intellectual disagreements which have provoked 
this work. 



APPENDIX A: INTEGRALS OVER 
DIELECTRIC FUNCTIONS 

Here we proof a very useful relation, which has 
been given in El. Therefore we consider the fol- 
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lowing integral including the dielectric function 



dm H(u>) 
dio / 1 



Ime~V) 
1 



4iri \uj + ir\ u> — ii] 



H(u)(f--f+) (Al) 



where /+ = 1 - 1/e and /~ = (/+)*. In the 
following we will assume that the function H{uS) 
is analytical and vanishes with ~ u>~ 2 for large u>. 
Since / ± (w) has no poles in the lower/upper half 
plane we have the identity 



^7^ = ^(0) (A2) 
2m (oj ± ir\) 



and all other combinations of with the denomi- 
nator vanish. If we would use the quantum dielec- 
tric function e we would have to add the residue of 
the poles at the Matsubara frequencies. Because 
we calculate only with the classical dielectric func- 
tion we can use (A.2) With t he h elp of the relation 
(A2) we compute easily for (Al) 



which proves relation (|70|). 



APPENDIX B: AN INTEGRAL 



Here we calculate the integral (174 



I[a, b] — dz 



l 

f i 

Z I . 2T 



dx- 



z 2 + 1 J a 2 x 2 + z 2 
o -1 



x I die z 1 cos(bzl 2 x) 



(Bl) 



The variable substitutions I — > p by p = ^fzl, x - 
z by z = yx and p — > e by P\fx = I leads to 

I[a,b] 

1 oo ^ ^ oo 

- 2 [dz I <///-r^ r / , [ dee~ y(? cos (be 2 ) 



dee ye cos (be 2 



j j y 2 x 2 + 1 a 2 + y 2 J 





V 2 / log(l + y 2 ) 



a 2 + y 2 



(B2) 



where the trivial x- integration has been carried 
out. The variable substitution e— >^ q by i/ye = 
q and y — > z by y = 1/z leads to 



I\«.b] = /dge-« 



OO 

y cos(6g 2 



l-z 2 log(l + 4 



z 2 + ^ 



(B3) 



Now we proceed an d us e an integral calculated in 
the next subsection 



dye^ 1 ^ os(1 + F) = 27r / dxx 2 <' 



1 ~2 
3 -A 



+a 7 re-/°(l+ 1 ° g(1 : a2) ) 



(B4) 



to obtain for (B3) 

7T3/2 1 + '"gC 1 "" 2 ) 



/[a, 6] 



4a 



7T 3 / 2 

2a 2 



da;- 



(^-x 2 )VT+te 



(B5) 



The last integrals is trivial and we end up with 



(A3) (0 



1. Another Integral 

Our task remains now to solve the integral 
l-y 2 log(l + ^) 



dye* 



y 



(B6) 



Because the complex function log(l + 1/y 2 ) has a 
cut from (0, i) we perform the integration along the 
path as depicted in figure [?] and write 



R 



-Cr + C r 



R r 

= 2iri Res 
= nae- c/a (l 



ii 



l-y 2 log(l + 4) " 

v*~T i ll,a 

y + a 2 

log(l-a 2 ). 



(B7) 



It is now easy to proof that in the limit r —> and 
i? — > oo the integration parts C vanish. Since the 
first two parts of (B7) represent just the desired 



integral I we have to calculate 

l-y 2 log(l + l ) 



= / dye tcy - 

I II i+r 



v 2 + ^ 
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\ lc , 1 -y 2l °g( 1 + F) + 27rz 
dye v 



y 



2 i J. 



-2tt / (te- 



rn 



Using (|B8|) and (|B7[ ) we obtain just (|B4|) 




FIG. 3. The complex integration path for the inte- 
gral |36l). 
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